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Abstract 
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^ \ A new time contour is used to derive real-time thermal field theory. Unlike previous 

•/^ ■ path integral approaches, no contributions to the generating functional are dropped. 

o\ : 
o . 
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In the path integral approach to equilibrium thermal field theories, one is calculating 
• thermal expectation values, typically 

> : 

'O ; Tc{ti,T2,...,Tn) = 

: Tr{e<^-PHTcMri)Mr2)...MrN)}/Tr{e^-PH}. (1) 

Here Tc indicates that the fields are path-ordered with respect to the order along a path, 
C, of their complex time arguments {r} |[ with appropriate sign changes when 
fermionic fields are involved . Throughout we shall suppress spinor and other indices as 
well as dependence on spatial variables where necessary. The results are not altered by 
their inclusion. The e^—f3H factor is used to shift the bra states in the trace by —i(3 
in time so that all times involved lie on a path, C, which must run in the complex time 
plane so that it ends —i(3 below its starting point. 
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Another limitation on C that is often imposed comes from studying the behaviour 
of the Green functions. Inserting complete sets of energy eigenstates into the thermal 
Wightman function, say Tai...ai^i where the fields are in a fixed order, we have 

rai...ajv(''"l5 •••5 "^n) 

= Tr{e^-(3Hct>aA'raMaAra,)-4aAra^)}lTr{e^-[3H}- 

= [eXp{-lEmAb-Ta^+Taj^)}. 

mi,m,2,...,m]sf 

(n$iV,=2ea;p{-zE„^.(-r„^, +r,^._J}) . ([] $iV,=i(m,#„^,(0)|m,+i))]/Tr{e<l>-/3iJ}. 

(2) 

The {\mj)} are a complete set of energy eigenstates for temperature 1 {tjii+n = ^i)- 
We assume that the energy sums in (Q) are uniformly convergent when the exponentials 
in (0) are of the form exp{zjEj} with ^e{zj} < 0. The thermal Wightman functions, 
raia2....a3({T}), are therefore bounded when the complex time arguments tj satisfy {r} G 
Aa^a2....a3 where 

- >Mv) >/m(r,^ (3) 

This then implies that the thermal Wightman function is also analytic in this region [Q. 

The path-ordered Green functions are merely a combination of products of these 
thermal Wightman functions (H) with path ordering theta functions. In order to ensure 
that the path ordered Green function rc({T}) is also analytic away from hypersurfaces 
of equal time arguments, one finds that one must ensure that the time curve C always 
has a decreasing imaginary part as one moves along C 0. This is a further condition 
that one imposes on the curve C (though a constant real part may be acceptable if one 
is prepared to accept generalised rather than analytic functions, see section 2.1 in 0). 

Given these two limitations on C one is then free to choose the time curve to suit the 
problem at hand. The standard choice is to choose one running parallel to the imaginary 
axis, C = C/ of fig.(0). This gives the ITF (Imaginary Time- Formalism) or Matsubara 
method. Because this has no real times present, one must make an analytic continuation 
to real times of the Green functions calculated if one wishes to obtain dynamical infor- 
mation. This operation is in principal highly non-trivial because the analytic structure 
at finite temperature is much more complicated than at zero temperature. 

To avoid this problem RTF (Real-Time Formalisms) were developed. The path in- 
tegral approach was highlighted by Niemi and Semenoff though RTF had also been 
considered earlier (see pf for a review). In these approaches one chooses the curve in the 
complex time plane to be C = Cr where Cr = Ci © C2 © C3 © C4 runs from — T + tq 
to +T + To (Ci), +T + To to +T + Tq — ia[3 (C2), +T + tq — ia[3 to — T + tq — ia[3 
(C3), and — T + tq — ia[3 to — T + tq — i[3 (C4) 0, |, |]. This is shown in Figure 0. The 
parameters a = 1 — a and tq is arbitrary and reflect some of the freedom of choice in the 
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path C. Physical results do not therefore depend on them. The limit of T ^ oo is taken 
to produce the standard RTF. 

The draw back of this approach is that one has to show that the vertical sections 
6*3,6*4 are irrelevant 0, H, |^. Otherwise one would have an unwieldy combination with 
few similarities to zero temperature Minkowskii space field theory, and one would not see 
the close relationship between this path integral method and the Thermo Field Dynamics 
and C*$*-algebra approaches to RTF (see for a review). To see that C*3,C*4 can be 
ignored, one has to take the T oo limit and then demand that the sources, j(r), are 
zero at 3?e{r} — cx) in manner reminiscent of the usual asymptotic conditions used in 
zero temperature field theory. 

Dropping these contributions is disquieting to many. The use of an asymptotic con- 
dition at infinite times in the past and future must be some sort of approximation in 
thermal field theories. Any real-time Green function is telling us about dynamics and 
therefore some sort of deviation from equilibrium. To use equilibrium theories in such 
problems, we are assuming that deviation from equilibrium is 'small' and can be ignored 
on 'short' time scales. After a finite time one must face the fact that the system is re- 
turning to equilibrium. This can not be done within equilibrium thermal theories. Thus 
one should not really be relying on what the theory is doing on infinite time scales. 

In practice, calculations of real-time Green functions in RTF and ITF seem to be 
in complete agreement when we compare the same type of Green functions (the same 
diagram usually corresponds to different Green functions in different formalisms 0, |^, 
The one area that still appears to be a problem is when one is interested in static (zero 
energy) Green functions such as appear in free energy calculations. In RTF, some of 
the relevant diagrams are infinite due to the appearance of singularities of the form 
[5{k^2 - m$2)]$A^ > 2 [|[ {pace g g]). Extra rules were shown to be needed for RTF 



calculations in these situations []10|, [TT[ . In the path integral derivation of these rules [IT 
the C*3,C*4 sections did contribute when static sources were involved. Only by using a 
simple trick leading to an extra Feynman rule could one use the standard RTF method. 

It would therefore be advantageous to find a path integral approach to RTF that did 
not involve dropping any parts yet gave the standard RTF Feynman rules. The starting 
point is to note that it is common to think of the C*i,C*2 of the standard RTF curve Cn of 
fig. (|1|) as having infinitesimal slopes downwards so that the path ordered Green functions 
defined on these curves are analytic. The difference in the imaginary part of the ends of 
the C*i,C*2 curves is also infinitesimal. This latter property is unnecessary for the Green 
functions to be well behaved and it is this condition which we drop. Thus we might try 
the new curve Cjy = Cm ®Cn2 shown in fig.(0). This curve satisfies all the requirements 
discussed above. It is convenient to parameterise the curve using a real variable, t, and 
a thermal label, a = 1, 2. We define a complex valued doublet, r$a, so that 

rMt) = {'Hl^"'^ ^^^^ = i (4) 
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The constants g^l, g^2 are 



Thus Cn runs from r$l(— aT) = —aT + I'yP + tq to T^l{aT) = aT + I'jaP/a (Civi) 
and then from r$2(a;T) = aT + Z'jaP/a to r<l>2(— aT) = —aT + i'~^(3 + tq {0^2)- The 
Tq is an arbitrary complex time and it is in the region of tq that we wish to investigate 
the Green functions. The ends of the curve must be separated by —i(3 so that we have 
7 + 7 = 1, though the Feynman rules will turn out to be completely independent of this 
7 parameter. We have also defined a + a = 1 which will turn out to play the same role 
for Cat as they did for Cr. 

The constant T is a time and one takes the limit T 00 so that the slope of the two 
sections Cni,Cn2 becomes infinitesimal as required. Also in this limit, Cn and Cr are 
then essentially identical in the region |r — ro| ^ T. The limit of T ^ 00 will therefore 
generate the usual RTF method and yet no part of the contour has been dropped. There 
can be no suspicion that RTF contains anything less than any other formalism such as 
ITF where a different curve is employed. 

It is enlightening to check this new approach to RTF by studying the simple example 
of a real relativistic self-interacting field theory with Lagrangian 

C = ^(9<l>/i0)<l>2 + im$20<l>2 - V[(f)]. 
The generating functional can be written as 

Z[j] = exp{V[-t-^]} . Zo[j] 

Zo[j] = exp{-^ dTdT^fj{T)A,{T - r$/)j(T$0} 

where 

- (D^ + m$2)Ac(r - r^f) = 4(r - r$/). (6) 

The Dc contains time derivatives taken along C. Also 6c{t — t^/) = 1 if r is further along 
C than r$/ otherwise 9c{t - r$/) = 0, ^6'c(r - r$/) = 6c{t - r$/) and 6c{t - r$/) = 1 
defines the delta function (c.f. [jT], |, |l2l). 

Using the KMS condition and the equal time commutation relation [0(t, x), (/)(r, y)] = 
as boundary conditions the solution of (||) is 

—t 1 

AJt, u) = — —\\e^-iuT+e^-6ue^+iuT]9JT) + \e^+iujT+e^-duje^-iuJT]9c(-T)} 

2uj 1 - e^-Pu 

(7) 

where u = |(A;$2 + m$2)$^| and we have taken the Fourier transform with respect to 
the spatial variables. 



Now consider C = Cn and define doublet fields and sources with real time arguments 
{t} in terms of the original fields and sources, whose time arguments lie on Cn, through 



a(t) = 0(r<l>a(t)); j$a(t) = ^<l>a j(r<l>a(t)). (8) 

The a = 1, 2 is the thermal label, and the r$a is the parameterisation of the new curve 
given in (^. The generating functional can then be written in the usual final RTF form 

Z[j] = expij: l_^^^aTdt g^aV[-^-^]}.Zo[J] 

Q, — 1,2 

Zo[j] = exp{-- V / ^aTdt dt^t i^a{t)/\^ah{t-t^t)i^h{t)} 

but where in this case no term has been dropped or absorbed into the normalisation. The 
two-by-two matrix propagator is defined to be 

A$a6(t - = Ac(r<l>a(t) - T^h{t^l)) 

The Feynman rules using these definitions lead to a factor of g<^a (as defined in (|^)) 
associated with each vertex. In the T — > cxd limit this leads to the usual additional 
factor of —1 appearing in type two vertices as compared with the type one vertices. 
The definition of the doublet fields and sources in (P) is fairly arbitrary and other useful 
definitions can be used which lead to variations in the Feynman rules [0]. 
When the limit T — oo is taken the doublet fields are then 

lim <P^l{t) = <P{t) 

lim 0$2(t) = 0(t - ial3) 

T — >oo 

where without loss of generalisation we have chosen ro=0. Thus we are indeed calculating 
the same Green functions as in standard RTF, for example 

d^NZlj] 

Tr{e<!>-l3HTMri)Mr2)...MrN)}/Tr{e^-l3H} 
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which is the connected time-ordered Green function of physical real-time fields. 

Calculating the Fourier transform of the propagator in this limit requires care with 
the various infinitesimal quantities. This is because 



lim / ^Tdt e$z(p + ie)t . 9(t)e^-iuj.(l - ir])t 

T— ►oo J-T 



= \l- \im{e^i{v-uj)Te^-{t + uori)T)]. (9) 

where p and u are real and e and t] are real and infinitesimal. Consider the calculation of 
A$ll. Then r/ = (1 — g^l)/i is positive and comes from the infinitesimal gradient given 
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to the otherwise horizontal Cni- As usual if only particle solutions are propagating at 
large times, this infinitesimal slope to the time path plays the same role as the usual le 
shift in the energy used to regulate the integral. However at non-zero temperatures we 
have both particle and anti particle solutions, and this means we also have terms where 
uj — *■ —u) in (|^). The slope of the time curve is of course the same so now we require 
that e > rj for a finite result. The same behaviour is found when looking at the particle 
contributions to the 9{—t) part of the propagator and so on for all parts of the matrix 
RTF propagator. In general we find that the infinitesimal slope given to time curves at 
non-zero temperature is not sufficient to guarantee convergence of the Fourier transform 
of the free propagator when T — oo. 

The problem is easily avoided. We are still free to give the time curves infinitesimal 
slopes but we must also have the usual le shift in the energies. We must then take the 
slopes to zero before we take the e — > limit. In terms of the RTF curves (both for this 
new approach and for the old style curves) we must take the T oo limit before the 
e — *• limit. 

With this proviso we find that the usual RTF propagators in momentum space are 
reproduced. In this example where the definitions have been chosen to give a factor of 
— 1 appearing in the type two vertices over the type one, we find that 

tA<l>ll(ko,uj) = + — 27T6(ko<l>2 - uj<^2) 

tA<^>12{ko,u}) = -^^^^^e^{a-^)(3ko/2 27r6{ko^2-uj^2) 

zA$21(fco,cu) = ^^^^e$-(a--)/3A;o/2 27r(5(A;o$2-cu$2) 
e^pu — 1 2 

tA^22{ko,to) = , + ^:^ -27r5(A;o$2 - ^$2) 

ko^2 - uj^2 - It e^fiuj - 1 

The physical origin of these problems is that T represents the time scale over which 
the system starts to come into equilibrium. In looking at dynamical information in the 
form of real-time Green functions, we have the inherent contradiction that we must be 
pushing the system out of equilibrium yet we assume the system is in equilibrium for 
all time. This is a good approximation on time scales less than T but is bad otherwise. 
Mathematical consistency therefore demands that this characteristic time scale is taken 
to infinity in such equilibrium schemes. It is not surprising that the T — > cxd limit must 
be taken first. 

This also explains why we should not be surprised if when using RTF to look at 
zero energy Green functions we have to work a little harder as these involve averages 
over all times including those of order T. This could be seen in the extra rule found 



to be necessary [10, 11|. When using the standard RTF time curve the vertical sections 



were shown to be giving all the contribution to the generating functional [TT[ but after 



some manipulations, it was shown the standard RTF approach could be used if one extra 
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Feynman rule was added. 

For the new time curve the formal manipulations of 0] work exactly as before. The 
ends of the curves are significant and so the infinitesimal slope can not be ignored. The 
same trick as before allows one to take account of this when using the new curve by 
introducing the same extra rule. Once again the new time curve completely reproduces 
the standard RTF method without dropping any contributions. This gives us further 
confidence that ITF and RTF contain the same physical information. 

I would like to thank R. Kobes, G. Kunstatter and the Institute for Theoretical Physics 
at the University of Winnipeg for their hospitality and the Central Research Fund of the 
University of London for financial support for my visit to Winnipeg. 
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Figure 1: The time paths used for RTF and ITF. 
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Figure 2 



The new time path for RTF. 
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